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Abstract We discuss time-dependent factorial cumulants in 
interacting nano-scale systems. Recent theoretical work has 
shown that the full counting statistics of non-interacting elec- 
trons in a two-terminal conductor is always generalized bi- 
nomial and the zeros of the generating function are conse- 
quently real and negative. However, as interactions are intro- 
duced in the transport, the zeros of the generating function 
may become complex. This has measurable consequences: 
With the zeros of the generating function moving away from 
the real-axis, the high-order factorial cumulants of the trans- 
port become oscillatory functions of time. Here we demon- 
strate this phenomenon on a model of charge transport through 
coherently coupled quantum dots attached to voltage-biased 
electrodes. Without interactions, the factorial cumulants are 
monotonic functions of the observation time. In contrast, as 
interactions are introduced, the factorial cumulants oscillate 
strongly as functions of time. We comment on possible mea- 
surements of oscillating factorial cumulants and outline sev- 
eral avenues for further investigations. 

Keywords Full counting statistics • noise • factorial 
cumulants • interactions • generalized master equations 

PACS 02.50.Ey • 72.70.H-m • 73.23.Hk 



1 Introduction 

The full counting statistics (FCS) of charge transfers in sub- 
micron electrical conductors has become an active field of 
research |[T]|2]|3]- Initially, investigations of FCS were pri- 
marily of theoretical interest, but several experiments EE 
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"ellTllFlSITSinnilinnilH have now clearly 

demonstrated that measurements of FCS are achievable and 
much progress has been made: Non-Gaussian voltage and 
current fluctuations have been measured in tunnel junctions ID 
l6l fT2l and quantum point contacts IfTSl . and the fourth and 
fifth current cumulants have been detected in an avalanche 
diode ifTSl . Additionally, real-time electron detection tech- 
niques ||5]|2l have paved the way for measurements of the 
FCS of charge transport in single ||9][l0][n][l6][I2][l8][l9] and 
double quantum dots ll8l fT4l . Following the initial measure- 
ments of the third cumulant of transport through quantum 
dots |l8]|9l, a series of experiments have addressed the condi- 
tional FCS ifTTI . the transient high-order cumulants III6III7I 
[T8l . and the finite-frequency FCS |fT9l in quantum dot sys- 
tems. The works on transient high-order cumulants showed 
that high-order cumulants generically oscillate as functions 
of basically any system parameter as well as the observation 
time lfT6l . 

Investigations of FCS are motivated by the expectation 
that more information about the fundamental transport mech- 
anisms can be extracted from the full statistical distribution 
of transferred charges than from the mean current and shot 
noise only |[l]|2]|3| . However, the fact that high-order cumu- 
lants generically oscillate makes it less clear exactly what 
information the high-order cumulants contain? In a recent 
work ||20| . we have been drawing attention to the use of fac- 
torial cumulants to characterize the FCS of charge trans- 
port in nano-scale electrical conductors. So far, factorial cu- 
mulants have only received limited attention in mesoscopic 
physics (but see Refs. fl21ll22ll23]l24ll ). However, as we have 
shown, the factorial cumulants never oscillate (unlike the or- 
dinary cumulants) for non-interacting two-terminal scatter- 
ing problems. This result is based on the recent finding that 
the FCS for non-interacting electrons in a two-terminal scat- 
tering setup is always generalized binomial and the zeros of 
the generating function for the FCS consequently are real 



2 



Dania Kambly, Christian Flindt 



and negative 12511261 : see Ref. IZTl for a discussion of multi- 
terminal conductors. In contrast, as interactions are intro- 
duced in the transport, the zeros of the generating function 
may become complex and the factorial cumulants start to 
oscillate f20l. This indicates that factorial cumulants may be 
useful to detect interactions among charges passing through 
a nano-sized electrical conductor 

The purpose of this work is to illustrate these ideas on 
a model of transport through coherently coupled quantum 
dots. In previous work [20;28], we considered systems de- 
scribed by classical master equations. We now turn to a sit- 
uation, where the quantum coherent coupling between two 
parts of the conductor is important. The system we consider 
is a double quantum dot (DQD) attached to external source 
and drain electrodes. We employ a generalized master equa- 
tion (GME) approach which allows us to treat strong cou- 
pling to the leads together with the coherent evolution of 
electrons inside the DQD. We treat two cases of particular 
interest: In the non-interacting regime, the DQD can accom- 
modate zero, one, or two electrons at a time, without addi- 
tional charging energy required for the second electron. We 
show that the factorial cumulants in this case do not oscillate 
as functions of the observation time and from the high-order 
factorial cumulants we extract the zeros of the generating 
function which are real and negative. Next, we consider the 
strongly interacting case, where double-occupation of the 
DQD is excluded. In this case, the time-dependent factorial 
cumulants oscillate - a clear signature of interactions in the 
transport - and the zeros of the generating function are com- 
plex. 

In the interacting case, we find that the Fano factor F, 
i. e. the ratio of the shot noise over the mean current, may 
either be super-Poissonian {F > 1) or sub-Poissonian {F < 
1). Super-Poissonian noise is typically taken as a signature 
of interactions in the transport [|2l, while no clear conclu- 
sion can be drawn from a sub-Poissonian Fano factor Inter- 
estingly, we find that the factorial cumulants may oscillate 
in both situations, showing that the factorial cumulants can 
provide a clear signature of interactions even when the cur- 
rent fluctuations are sub-Poissonian. We conclude our the- 
oretical investigations of factorial cumulants by examining 
the influence of dephasing of electrons passing through the 
DQD II29II30L for instance due to a nearby charge detector. 

The paper is organized as follows: in Sec. |2]we intro- 
duce the essential terminology used in FCS and provide the 
basic definitions with a special emphasis on factorial cu- 
mulants and the concept of generalized binomial statistics. 
In Sec. [3] we then turn to the asymptotic behavior of high- 
order cumulants (both ordinary and factorial cumulants) and 
show why the high-order factorial cumulants do not oscillate 
for transport of non-interacting electrons in a two-terminal 
conductor. In Sec.|4]we introduce a model of electron trans- 
port through a DQD described by a Markovian GME, while 



Sec.|5]is devoted to the details of our calculations of time- 
dependent factorial cumulants. In Sec. |6] we demonstrate 
how interactions on the DQD give rise to clear oscillations of 
the high-order factorial cumulants with the zeros of the gen- 
erating function moving away from the negative real-axis 
and into the complex plane. Finally, Section Q is dedicated 
to a summary of the work as well as our concluding remarks. 

2 Full counting statistics & factorial cumulants 

Full counting statistics concerns the quantum statistical pro- 
cess of electron transport in mesoscopic conductors. The 
full counting statistics (FCS) is the probability P{n,t) that n 
electrons have traversed a conductor during a time span [Q,t] 
of duration t. The information contained in the probability 
distribution may equally well be encoded in the generating 
function (OF) defined as 

^{z,t)='£Pin,t)z". (1) 

The normalization condition for the probabilities, ^„ P{n,t) = 
1, implies for the GF that ^^(z = l,f) = 1. Important infor- 
mation about the charge transport can be obtained from the 
GF: If the transport process consists of several independent 
sub-processes, the GF factors into a product of the GFs cor- 
responding to each of these sub-processes, similarly to how 
the partition function in statistical mechanics may be writ- 
ten as a product of the partition functions for each indepen- 
dent sub-system. Elementary transport processes can thus be 
identified by factorizing the GF. In the case of transport of 
noninteracting electrons through a two-terminal conductor, 
Abanov and Ivanov have shown recently that the GF can be 
factorized into single-particle events of binomial form ll25l 
l26ll . Such distributions have been dubbed generalized bino- 
mial statistics Il31ll25ll26l . which will be of central impor- 
tance in this work. 

Several useful statistical functions and quantities can be 
obtained from the GF. First, we can define a moment gener- 
ating function (MGF) 

^(z,0 -^(e-',f), (2) 

which generates the statistical moments of n by differentia- 
tion with respect to the counting field z at z = 

(«"')(f) = 01.^0 = i:«"T(n,f). (3) 

The MGF of a transport process composed of several inde- 
pendent processes factors into a product of the correspond- 
ing MGFs. However, the moments of the full process are 
not related to the moments of the individual sub-processes 
in a simple way. This motivates the definition of cumulants. 
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also known as irreducible moments. The cumulant generat- 
ing function (CGF) is defined as the logarithm of the MGF 



y{z,t)=\og[^{z,t)]=\og[^{e\t)], 



(4) 



which again delivers the cumulants of n by differentiation 
with respect to z at z = 0: 



((«'"))(r) = af^(z,f)|,^o. 



(5) 



The first cumulant is the mean of n, ((«)) = (n), the second 
cumulant is the variance, ((n)) = (n^) — (n)", and the third is 
the skewness, ((n^)) = ((«— («))^). Itis easy to show that the 
cumulants of a transport process are simply the sum of the 
cumulants corresponding to each independent sub-process. 
Moreover, for a Gauss distribution only the first and second 
cumulants are non-zero, while all higher cumulants vanish. 
In this respect, one may use cumulants of a distribution as a 
measure of (non-)gaussianity. 

The conventional moments and cumulants, as defined 
above, have been investigated intensively in the field of FCS 
Il3l- The zero-frequency cumulants of the current are given 
by the long-time limit of the cumulants of n as 



((«'"» 



(6) 



The current is treated as a continuous variable and continu- 
ous variables are typically characterized by their cumulants. 
However, another interesting class of statistical quantities 
exists, which has received much less attention in FCS. These 
are the factorial moments and the factorial cumulants, which 
are mostly discussed in the context of discrete variables 1321 
[33)1 . The number of counted electrons « is obviously a dis- 
crete variable, and it is natural to ask if the current cumu- 
lants, as defined in Eq. (|6]l, carry signatures of this discrete- 
ness? 

The factorial moments are again generated by a factorial 
MGF, which can be defined based on the GF in Eq. ([T]i- The 
factorial MGF is defined as 

^f(z,f) =^(2+1,0 (7) 

and the corresponding factorial moments read 

(«"') F [t) = d'^J^F {z,t)\,^Q = («(n - 1 ) • • • (« - m + 1 )) (8) 

in terms of the ordinary moments. In analogy with the con- 
ventional CGF, the factorial CGF is defined as 

^f(z,f)=logKf(z,f)] =log[^(z+l,f)] (9) 

and the corresponding factorial cumulants read 

(r) = d'^^F {z, f ) Iz^o = ((«(« - 1 ) • • • (« - m + 1 )» . 

(10) 



As mentioned above, factorial moments and factorial cumu- 
lants are of particular interest when considering probability 
distributions of discrete variables. For example, for a Pois- 
son process with rate F, which is the physical limit of rare 
events, the FCS is well-known and reads 



n ' 



(11) 



The corresponding GF then becomes 

^(z,f) ^e^'*^'"'', (Poissonprocess) (12) 

from which it is easy to show that the cumulants are 

((«'"))(?) =rf, (Poissonprocess) (13) 

for all m ~ 1,2, . . .. In contrast, the first factorial cumulant 
reads 



((«))/r(f) = Ff , (Poissonprocess), 

while all higher factorial cumulants are zero 

((«'")) /r(f) = 0, m > 1 (Poissonprocess). 



(14) 



(15) 



Thus, similarly to how ordinary cumulants are useful as mea- 
sures of (non-)gaussianity, we may use factorial cumulants 
to characterize deviations of a distribution from Poisson statis- 
tics. It is also clear that a factorial cumulant of a given order 
is the sum of the factorial cumulants of all independent sub- 
processes, in the same way as for the cumulants. 

Throughout this work we will rely on an important re- 
sult by Abanov and Ivanov, who have shown that the FCS of 
non-interacting electrons in a two-terminal scattering prob- 
lem is always generalized binomial II25II26I . In this case, the 
GF takes the special form 



ryy/ \ generalized _ 

^{z,t) = Z 

binomial 



n(i 



-Piz), 



(16) 



where the (time-dependent) /7,'s are real with < p, < 1- 
Each factor in the product can be interpreted as a single bi- 
nomial charge transfer event occurring with probability p,. 
The factor in front, z^^, corresponds to a deterministic back- 
ground charge transfer 



(17) 



opposite to the positive direction of the mean current. For 
uni-directional transport, ^, pj = («) and 2 = 0, whereas Q 
is non-zero for bi-directional transport due to thermal fluctu- 
ations for example. For uni-directional transport, the (time- 
dependent) Fano factor then reads 



F{t) = generalized_ L'PKl ' Pi) 



(«)(f) binomial 



Li Pi 



(18) 



which is always smaller than unity, corresponding to a Pois- 
son process. Following this reasoning, a super-poissionian 
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Fano factor, F > I, can be taken as a sign of interactions. 
Super-Poissonian noise was recently measured in an exper- 
iment on transport of interacting electrons through a double 
quantum dot lf34l . Still, the noise may also be sub-Poissonian, 
F < 1, in the presence of interactions. 



3 High (factorial) cumulants 

In this Section we are interested in the generic behavior of 
high-order (factorial) cumulants. To this end, we first discuss 
an approximation of high derivatives [35, 36] and then show 
how these ideas can be applied in the context of FCS. 

In the following, we consider a generic (factorial) CGF 
'$^{F) ^nd assume that it has a number of singularities zj in 
the complex plane. Close to each of these singularities, we 
may approximate the (factorial) CGF as 



for z: 



(19) 



for some constants Aj and jij. Here, the constant jij is deter- 
mined by the nature of the singularity, for instance pij~~^/2 
corresponds to a square-root branch point, while an integer 
value of pLj would correspond to the order of a pole at z = Zj- 
Using the first Darboux approximation li35ii36jlT6l . we may 
now evaluate the (factorial) cumulant of order m by differ- 
entiating the expression in Eq. ( fT9] l m times with respect to z 
at z = and sum over the contributions from all singularities 
as 



AjB„ 



{F) 



(20) 



Here we have introduced the polar notation zj = |z;|e'"^''^ 
together with the factors 



B„ 



■V-j 



(21) 



Equation ( l20l i is particularly useful if the sum can be re- 
duced to only a few terms. For high orders, the singularities 
closest to z = dominate the sum, which leads to a consider- 
able simplification. For example, if a single complex conju- 
gate pair of singularities, zo = Izol'''''^^"''! andzo = IzoT""'^'-'''!, 
are closest to z = 0, the high-order (factorial) cumulants can 
be approximated as 

((«"'» (n ^ ^l^iSr '^"^ [("^ + ^--g^o - ^g^o] . (22) 

This result shows that the absolute value of the (factorial) cu- 
mulants generically grows factorially with the cumulant or- 
der m, due to the factors B,,,,^,,, and that they tend to oscillate 
as a function of any parameter, including time f , that changes 
argzo- Such universal oscillations have been observed exper- 
imentally in electron transport through a quantum dot |1T6] 

mm- 



In contrast, in the particular situation, where there is just 
a single dominant singularity zo on the real-axis, the high- 
order (factorial) cumulants can be approximated as 



((«'"»( 



"I. Mo 



1^0 



|"i+Mo ■ 



(23) 



In this case, the factorial growth with the order persists, but 
no oscillations are expected as long as the dominant singu- 
larity Zo stays on the real-axis. 

Let us now consider non-interacting electrons in a two- 
terminal conductor As shown by Abanov and Ivanov ll25l 
26], the statistics is generalized binomial in this situation 
and the GF takes on the form given by Eq. (fTST l. The corre- 
sponding cumulants are complicated functions of the prob- 
abilities Pi. In contrast, the factorial cumulants are simply 



^^^^,„^^^ generalized (_ ^ 1 ^ _ ^ ) , 
binomial 



(24) 



For uni-directional transport {Q = 0), the largest probability 
/'max will dominate the high factorial cumulants, which can 
be approximated as 

^(-l)'"-'(m-l)!Kax- (25) 

This expression can also be understood from Eq. (l20t by 
noting that the factorial CGF has logarithmic singularities 
at values of the counting field z, where the factorial MGF 
is zero. Combining Eqs. (|7]i and (fTSI l. we easily see that the 
factorial MGF has zeros at Z; = — 1 / Pj < — 1 . Moreover, the 
zero corresponding to the largest probability pmax is closest 
to z = and will dominate the high factorial cumulants as 
seen in Eq. (l25T l. 

We have seen above that high-order cumulants tend to 
oscillate as functions of basically any parameter, with or 
without interactions. In contrast, as our analysis also shows, 
the high-order factorial cumulants never oscillate for non- 
interacting electrons in a two-terminal scattering problem. 
This behavior can be traced back to the factorization of the 
GF in Eq. (fTSI l. which implies that the singularities of the 
factorial CGF are always real and negative. This makes fac- 
torial cumulants promising candidates for the detection of 
interactions in FCS. In particular, oscillating factorial cumu- 
lants must be due to interactions. In our previous work iBOl . 
we employed these ideas to incoherent electron transport 
through a single quantum dot. We showed how interactions 
may cause the singularities of the factorial CGF to move 
away from the real-axis and into the complex plane, making 
the high-order factorial cumulants oscillate. 

In the following we apply these ideas to electron trans- 
port through a DQD, where the electrons may oscillate quan- 
tum coherently between the two quantum dots. In contrast 
to our previous work lEOl . we consider not only the time- 
dependent factorial cumulants of the transferred charge, but 
also the zero-frequency factorial cumulants of the current. 
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4 Coulomb blockade quantum dots 

We consider two-terminal nano-scale conductors connected 
to source and drain electrodes. Charge transport is described 
using a generalized master equation (GME) for the reduced 
density matrix p of the conductor obtained by tracing out the 
electronic leads. The GME accounts for the coherent evolu- 
tion of charges inside the conductor as well as the transfer 
of electrons between the conductor and the leads. To evalu- 
ate the PCS, it is convenient to unravel the reduced density 
matrix with respect to the number of electrons n that have 
been collected in the drain electrode during the time span 
[0,f] 07II38I . With this n-resolved density matrix p{n,t) at 
hand, the PCS is obtained by tracing over the states of the 
conductor 



P(«,r)=Tr[p(«,f)]. 



(26) 



Similarly, we recover the original reduced density matrix by 
summing over n, i. e. 



p(0 = Lp(«,0- 

n 

From these definitions, the GF reads 
$#(z,0=i:Tr[p(«,0]z"=Tr[p(z,f)], 



(27) 



where we have introduced the z-dependent reduced density 
matrix 



(29) 



The particular conductor we now discuss consists of two 
quantum dots in series attached to source and drain elec- 
trodes. A schematic of the system is shown in Pig. [T] The 
inter-dot Coulomb interaction can be tuned, so that both 
regimes of noninteracting and interacting electrons can be 
reaUzed and compared. Disregarding the electronic spin de- 
gree of freedom (for example due to a strong magnetic field), 
the double quantum dot can be occupied by either zero, one, 
or two (additional) electrons. Experimentally, the charge on 
the quantum dots can be measured using a nearby quan- 
tum point contact (QPC), whose conductance it sensitive to 
the occupations of the individual quantum dots 191 1141 . This 
charge detection protocol makes it possible to deduce the 
number of electrons that have passed through the DQD in 
a given time span. If the QPC is not sensitive to the charge 
occupations of the individual quantum dots, but only to the 
total charge, the measurement is not expected to destroy 
the coherent oscillations between the quantum dots. On the 
other hand, if the QPC measures the charge states of the in- 
dividual quantum dots, it introduces decoherence in the dy- 
namics of electrons inside the DQD 13911401 . 

The full many-body Hamiltonian of our system reads 



H — HdqD + ^^leads + Hj + Hqpc + HoQD-QPC- 



(30) 




QPC 



Fig. 1 Double quantum dot (DQD) coupled to a QPC charge detec- 
tor The upper panel shows the DQD. The tunnel coupling between the 
quantum dots is denoted as £2 and Ei and En are the single-particle 
levels of the two quantum dots with detuning e = Er — E^. The tun- 
neling rate from (to) the right (left) quantum dot to (from) the right 
(left) reservoir dot is denoted as -rj((L). The lower panel shows the QPC 
charge detector, which measures the charge occupation of the DQD. 
The QPC may couple asymmetrically to the DQD such that the charge 
occupation of the individual quantum dots can be resolved. 



(28) consists of the Hamiltonian of the DQD 



Hdqd = EidldL + EpdgdR + Q {d^dp + ajjOz,) + U fiLnp, 

where the operators fl| and d'j^ create an electron on the left 
and right quantum dot level with energy El or Ep, respec- 
tively. The tunnel coupling between the levels is denoted as 
£2 and = «««« = 0,1, a ~ L,R is the occupation number 
operator of each quantum dot. The inter-dot Coulomb inter- 
action is denoted as U. The electrons in the leads are treated 
as non-interacting and are given by the Hamiltonian 



^leads 



k.a=L,P 



(31) 



where the operators fl|^ create an electron in lead a = L,R 
with momentum k and energy e^a . The coupling between the 
DQD and the leads is accounted for by the standard Hamil- 
tonian 



Ht= (fkadlo,da + tladlcdka), 



(32) 



k.a=L,P 



which connects the left (right) lead to the left (right) quan- 
tum dot. Pinally, the QPC is modeled as a tunnel barrier 



H, 



QPC 



ka ka 



Y^{hk'clLCi,in + tl^,clj^Ci^,l), 



k,a=L.R 



k.k' 



where the first sum corresponds to the electronic reservoirs 
on the left {a = L) and right side (a = /?) of the QPC and 
the second sum describes the coupling of states in different 
leads with tunnel coupUng t^jft. 
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If the QPC only couples to the total charge of the DQD, 
the charge detection is not expected to cause decoherence of 
the coherent oscillations of electrons inside the DQD. It is, 
however, interesting to investigate how an asymmetrically 
coupled QPC will affect the transport in the DQD. To this 
end, we assume that the QPC, besides the coupling to the 
total charge, has an additional capacitive coupling to the left 
quantum dot only. The charge occupation of the left quan- 
tum dot modulates the transparency of the QPC according 
to the Hamiltonian 

HdQD-QPC = «i(^^H'cii4fl + ^Kk'^lR^k'l)- (33) 

k,k' 

Here dtkk' is the change of the QPC tunnel coupling in re- 
sponse to an (additional) electron occupying the left quan- 
tum dot. 

We now follow Gurvitz in deriving a Markovian GME 
for the reduced density matrix Pdqd of the DQD obtained by 
tracing out the electronic leads and the QPC. The details of 
the derivation can be found in Refs. P9II40I . We assume that 
a large bias is applied between the electronic leads, such that 
electron transport is uni-directional from the left to the right 
electrode. The electronic reservoirs have a continuous den- 
sity of states and the discrete levels of the DQD are situated 
well inside the transport window. Under these assumptions, 
we may formulate a Markovian GME for the «-resolved re- 
duced density matrix PDQD{n,t). Its diagonal elements are 
the probabilities for the DQD to be either empty, having 
only left or right quantum dot occupied, or to be doubly- 
occupied, while « electrons have been collected in the right 
lead during the measuring time t. 

The diagonal elements of PDQD{n,t) are denoted as po{n,t) 
piin^t), pi{{n,t), and Pci{n,t). Additionally, we need the co- 
herences between the left and the right quantum dot lev- 
els, denoted as PLfl(n,0 and PRL{n,t)- Coherences between 
states with different charge occupations are excluded. Since 
the off-diagonal elements fulfil PRL{n,t) = plj^{n,t), it suf- 
fices to consider the real and imaginary parts of PLR{n,t). 
The elements of the reduced density matrix can then be col- 
lected in the vector 



|p(«,f))) = \pO,PL,PR,'S{[pLR],Z[pLR],PdV{fht). 



(34) 



The corresponding z-dependent reduced density matrix fol- 
lows from the definition in Eq. ( |29] l and reads 



ip(z,o))^i:ip(«,o))^"- 

n 

The Markovian GME then takes the form 
a,|p(z,f)))=M(z)|p(z,0)), 



(35) 



with the rate matrix reading 



Mil) 



-a 


zFr 











a 








-2X2 


zFr 





-Fr-Tl 





2£2 











r 


-e 





Q 


£2 


e 


r 








a 








-Fr 



(37) 



where e = Er — Ei is the energy detuning of the two quan- 
tum dot levels. Additionally, the rate 



r = -{rR+a + r) 



(38) 



determines the decay of the off-diagonal elements of Pdqd («, 
and the broadening of the electronic levels. The electronic 
tunneling rates depend on the charge occupation of the DQD 
and read 



Pa 



and 



27Z 

— &a{Ea)\tka\ 



2n 



a 



L,R 



ra^=i^^a{Ea + U)\tkoc\^, a ^ L,R, 
11 



(39) 



(40) 



where denotes the density of states in lead a=L,R, and 
the tunneling amplitudes tka are assumed to be ^-independent. 
Here, II is the tunneling rate from the left lead onto the left 
quantum dot, if the DQD is empty initially. On the other 
hand, if the right quantum dot is already occupied, electrons 
tunnel into the left quantum dot at rate T/,. Similarly, elec- 
trons tunnel from the right quantum dot into the right lead 
with rate Fr, if the left quantum dot is empty, and with rate 
Fr, if the left quantum dot is occupied. Without inter-dot in- 
teractions, U = 0, we have Il(^) Pl{r)- Factors of z have 
been included in the off-diagonal elements of M(z) corre- 
sponding to charge transfers from the right quantum dot to 
the right lead. 

Finally, the decoherence rate introduced by the QPC is 
given by |401 



(41) 



where Vd is the bias applied accross the QPC. The transmis- 
sion probability for electrons to tunnel through the QPC is 



r = (2;r)2|4_,,|2^i^«, 



(42) 



when the left quantum dot is empty. In contrast, when the 
left quantum dot is occupied, the QPC transmission reads 



i27zf\tj,k + 8ij,k\^^L^R. 



(43) 



(36) 



Above, the symbols !^l(r) denote the density of states in the 
left (right) lead of the QPC. 
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Fig. 2 Fano factor versus the energy detuning e. Results are shown 
with (U ^ 0) and without (U = 0) strong Coulomb interactions on the 
DQD. With strong Coulomb interactions, only or 1 electrons can 
occupy the DQD. In contrast, without Coulomb interactions the DQD 
may also be doubly occupied. The QPC is coupled symmetrically to 
the DQD (7 = 0). The other parameters are fj; = i^Fi and Q = hFi. 
The squares and stars mark e = 0.812 and e = 2.2Q, respectively. 
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Fig. 3 Time-dependent factorial cumulants without interactions (U = 
0). The factorial cumulants ((;i"'))(f) of order ;;j = 11 through m = 14 
are shown as functions of time. The results correspond to the point 
marked with a red star in Fig.|2] The factorial cumulants do not oscil- 
late, as expected without interactions. The full lines indicate numerical 
results, while circles show the approximation given by Eq. i25\ . 



5 Calculations 

We now evaluate the FCS by formally solving Eq. (|36] |. 
We consider fluctuations in the stationary state, which we 
suppose has been reached at f = 0, when we start count- 
ing charges. The stationary state is denoted as |0)) and is 
obtained by solving M(z = 1)|0)) =0 with the normaliza- 
tion condition ((0|0)) = 1, where ((0| = [1, 1, 1,0,0, 1]. From 
Eq. ( |36] |, the GF can now be written as 

^^(z,r)-((OkM«'|0)). (44) 

It is easy to verify that this expression fulfils the condition 
= l,f) = 1 for the GF. It is a general, formally exact 
result, which yields the complete FCS at any time given 
the matrix M(z). However, in practice the expression may 
be difficult to evaluate due to the matrix exponentiation, in 
particular if the aim is to calculate the high (factorial) mo- 
ments or (factorial) cumulants. In our recent work lEOl . we 
developed a simple method to evaluate the high-order, time- 
dependent statistics for these types of problems and we will 
also be using this method here. For details of the method, 
we refer the interested reader to Appendix A of Ref . [l20| . 

In addition to the finite-time FCS, it is interesting to in- 
vestigate the GF at long times. In this limit, the GF takes on 
a large-deviation form 

^(z,f)oce'®W, (45) 

where the rate of change is determined by the eigenvalue of 
M(z) with the largest real-part, i. e. 

&{z)=max{Xj{z)}. (46) 
J 



From (z) we may obtain the (factorial) CGF for the zero- 
frequency cumulants of the currents. The zero-frequency cu- 
mulants of the current are 

((/"')) =5-0(^-^)1,^0, (47) 
while the corresponding factorial cumulants read 

((r))f = 5f0(z+ 1)1,^0. (48) 

In general, we can assume that the matrix M(z) at z = 1 
has a single eigenvalue equal to zero, i. e. Ao(l) = 0, cor- 
responding to the (unique) stationary state, while all other 
eigenvalues have negative real-parts, ensuring relaxation to- 
ward the stationary state. For values of z close to unity, we 
expect that Ao(z) develops adiabatically from and still de- 
termines @{z) such that 0(z) = Ao(z) forz — 1. The deriva- 
tives of Ao(z) with respect to z at z = 1 then determines the 
(factorial) cumulants of the current according to Eqs. (|47| 
andgH). 

Again, for large matrices M(z), it might not be viable 
to directly calculate the eigenvalue Ao(z) and its derivatives 
with respect to z at z = 1 . This problem may be circumvented 
by considering the calculation of Ao(z) as a perturbation 
problem around z = 1 . The matrix M(z) is written as M(z) = 
M(l) + 5M(z), where M(l) is the unperturbed matrix with 
eigenvalue Ao(z = 1) = and 5M(z) = M(z) - M(l) is the 
perturbation. The eigenvalue Ao(z) can then be calculated 
order by order in z using the recursive perturbation method 
developed in Refs. 1411142114311 . This method yields the (fac- 
torial) cumulants of the current and will be used below. 

Finally, it is important to understand the connection be- 
tween the FCS at finite times and in the long-time limit. As 
discussed in the previous section, the high (factorial) cumu- 
lants are related to the singularities of the (factorial) CGF in 
the complex plane of z- At finite times, the (factorial) CGF 
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Fig. 4 Time-dependent factorial cumulants with strong Coulomb interactions. The factorial cumulants {{n'"))f(t) of order m = 11 through m = 14 
are shown as functions of time. The panels correspond to the points marked with a blue square (a) and a blue star (b) in Fig. |2] The curves are 
shifted for clarity. We show the logarithm of the absolute value of the factorial cumulants such that downwards-pointing spikes correspond to a 
factorial cumulant going through zero. Full lines are numerical results, whereas empty circles correspond to the approximation given by Eq. j22b . 



has logarithmic singularities at values of z, where the (facto- 
rial) MGF is zero. In contrast, in the long-time limit, the sin- 
gularities of the (factorial) CGF are determined by the singu- 
larities of the eigenvalues of M(z). Typically, the eigenval- 
ues have square-root branch points at the degeneracy points, 
where two eigenvalues are equal, i. e. Ao(zc) = ^\{zc) for 
some Zc- Considering now the GF at finite times close to 
such a degeneracy point, we may approximate the GF in 
Eq. (lUli as 



-ci(z> 



(49) 



where the coefficients Cj (z) depend on the initial condition 
and only the contributions from the two largest eigenvalues 
have been included. Solving for the zeros of we ob- 

tain the equations 



Ao(z)=Ai(z) + 



log{ci(z)/co(z)}-f /;r(2«-f-l) 



(50) 



where n is an integer Importantly, we see that the second 
term on the right hand side vanishes in the hmit t ^ oo. This 
analysis shows that the zeros of the GF as functions of time 
move towards the solutions of the equation h){z) = Ai(z), 
which also determines the branch-point singularities in the 
long-time limit cf. the discussion above. This connects the 
finite-time FCS with its long-time behavior 



6 Results 

We are now ready to illustrate the use of factorial cumu- 
lants on the concrete example of charge transport through 
a DQD. We analyze several different parameter regimes of 
the system which are discussed in turn. To begin with, it is 



instructive to consider the Fano factor F of the transport in 
the long-time limit 



F = 



(n){t) 



(I) 



(51) 



given as the ratio of the zero-frequency current noise over 
the mean current. Figure |2] shows the Fano factor as a func- 
tion of the energy dealignment e without any decoherence 
due to the QPC, 7 = 0. We present results with (U ^ 0) and 
without interactions {U = 0). For the non-interacting case, 
the Fano factor is never super-Poissonian (F > 1) as ex- 
pected for uni-directional transport with generalized bino- 
mial statistics. In contrast, for the interacting case there are 
certain ranges of the dealignment, where the noise becomes 
super-Poissonian. However, there is also a range of dealign- 
ments around e = 0, where the noise is sub-Poissonian (F < 
1), and in this regime a measurement of the Fano factor 
would not give any clear evidence of interactions. We note 
that recent noise measurements on transport through ver- 
tically coupled quantum dots are in qualitative agreement 
with the results shown for the interacting case Il44j|34ll45l . 

We mark two points on the curves corresponding to val- 
ues of the dealignment, where the noise in the interacting 
case is either sub-Poissonian (squares) or super-Poissonian 
(stars). As we demonstrate now, the factorial cumulants give 
clear signatures of the interactions even in the cases with 
sub-Poissonian noise, where no conclusions can be drawn 
from the Fano factor alone. 

In Fig. [3] we show the time-dependent factorial cumu- 
lants for the non-interacting case, corresponding to the point 
marked with a star in Fig. |2] For the point marked with a 
square similar results are obtained. Without interactions, the 
FCS is generalized binomial and the factorial cumulants are 
expected to follow Eqs. ( l24b and ( |25] ), which predict no os- 
cillations of the factorial cumulants as functions of time or 
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Fig. 5 Motion of the dominant singularities with time. The singular- 
ities have been extracted from the high order factorial cumulants in 
Figs. [3] |4] and [6] The blue stars and squares are complex singularities 
corresponding to the oscillations of the factorial cumulants in Fig. |4] 
The red stars and squares are real singulaiities corresponding to the 
non-interacting case considered e. g. in Fig. |3] The encircled points 
correspond to the long-time limits considered in Fig. |6] The time is 
varied from = li/Q to tf = 50h/Q in steps of At = h/Q. 



any other parameter. This prediction is confirmed by Fig. |3] 
where a clearly monotonic behavior is found as a function 
of time. Moreover, from the calculated factorial cumulants, 
we may extract pmax in Eq. ( l25l l as a function of time. Insert- 
ing pmax back into Eq. dZSl ). we can compare this expression 
with the numerical results for the high-order factorial cumu- 
lants. In Fig. [3] the predicted behavior based on Eq. (l25T l is 
shown with circles and is seen to be in very good agreement 
with the full numerics. 

Next, we turn to the time-dependent factorial cumulants 
in the interacting case. In Fig. |4] we show the high-order 
factorial cumulants corresponding to the two dealigments 
marked with stars and squares in Fig. |2] In this case, the 
factorial cumulants oscillate as functions of time in con- 
trast to the non-interacting situation, where no oscillations 
are observed. To best visualize the oscillations, we show the 
logarithm of the absolute value of the factorial cumulants. 
Downwards-pointing spikes then correspond to the factorial 
cumulants passing through zero and changing sign. The os- 
cillating factorial cumulants are a clear signature of inter- 
actions in the transport and they show that the FCS for this 
system is not generalized binomial, neither when the noise 
is sub-poissionian (square) nor super-poissionian (star). 

Again, we can understand the high-order factorial cumu- 
lants using the expressions from Sec. |3] In this case, when 
the FCS is not generalized binomial, the high-order facto- 



rial cumulants are expected to follow Eq. (l22t . which as- 
sumes that the factorial CGF has a complex-conjugate pair 
of singularities. At finite times, the factorial FCS has loga- 
rithmic singularities corresponding to the zeros of the facto- 
rial MGF With only a single dominant pair of singularities, 
Zq and Zq, the expression for the high-order factorial cumu- 
lants, Eq. ( |22] | simplifies for finite times to III6II20I 



((«'")) 



2{m-l) 

i 

kol 



■cos (marg [zq] 



(52) 



From four consecutive high-order factorial cumulants, we 
can solve this relation for the dominant pair of singular- 
ities as functions of time using the methods described in 
Refs. l46B3]|20ll (see e. g. Appendix B of Ref. O). Insert- 
ing the solution back into Eq. ( l52b , we can benchmark the 
approximation against the numerical results. The approxi- 
mation is shown with circles in Fig. |4] and is seen to be in 
excellent agreement with the full numerics. 

Having extracted the dominant singularities as functions 
of time in the non-interacting and interacting cases, we may 
investigate their motion in the complex plane. In Fig. |5]we 
show the dominant singularities as they move with time. In 
the non-interacting case {U = 0), corresponding to the fac- 
torial cumulants shown in Fig. |3] the dominant singularity 
(marked with a red star) moves along the negative real-axis 
as expected for generalized binomial statistics. This behav- 
ior should be contrasted with the interacting case (U > 0), 
corresponding to the factorial cumulants shown in Fig. |4] 
In this case, the dominant singularities (marked with blue 
squares and stars) are no longer real and they now move 
in the complex plane as functions of time. We stress that 
this behavior cannot occur for a non-interacting system and 
should thus be taken as a signature of interactions. 

In Fig. |5] we also indicate the points in the complex 
plane to which the dominant singularities move in the long- 
time limit (encircled points). As discussed in the previous 
section, these points correspond to the dominant singulari- 
ties of the factorial CGF for the zero-frequency factorial cu- 
mulants of the current. We extract the position of these sin- 
gularities by calculating the high order factorial cumulants 
of the current ((/'"))f using the recursive scheme developed 
in Ref. here adapted to the calculation of factorial cu- 
mulants. The results for the factorial cumulants of the cur- 
rent are shown in Fig.|6]as functions of the order m. Together 
with the numerical results, we show the approximation in 
Eq. ( l22b with full lines. From four consecutive high-order 
factorial cumulants we have extracted the parameters enter- 
ing Eq. (I22I) using the method proposed in Ref. 1431 . Typi- 
cally, in the long-time limit, the singularities are square-root 
branch points such that = —1/2 in Eq. (l22l i. Figure |6] 
shows that Eq. (l22l i provides an excellent approximation of 
the numerical results and it allows us to extract the dom- 
inant singularities in the long-time limit (encircled points) 



10 



Dania Kambly, Christian Flindt 




5 10 15 20 25 

m 



Fig. 6 Factorial cumulants {{I"'))f of the current as functions of their 
order m. We coinpare numerical results (inarked with symbols) and the 
approximation (full line) given by Eq. \22\ . Panel (a) coiresponds to 
the point marked with a blue square in Fig. |2] Panel (b) corresponds to 
the point marked with a blue star. The corresponding dominant singu- 
laiities extracted from the numerical results are encircled in Fig.|5] 



in Fig.|5] As anticipated, the dominant singularities at finite 
times move towards the long-time singularities. In the long- 
time limit, the singularities cease to move with time and the 
high-order (factorial) cumulants will no longer oscillate as 
function of time (but still as functions of other parameters). 

Finally, we turn to the influence of detector-induced de- 
phasing. In Fig.|2]we consider the situation where the QPC 
charge detector is asymmetrically coupled to the two quan- 
tum dots, thereby causing dephasing of electrons passing 
through the DQD. Due to strong Coulomb interactions, the 
DQD can only be either empty or occupied by one (addi- 
tional) electron at a time. Without detector-dephasing (7 = 
0), clear oscillations of the tenth factorial cumulant of the 
current are observed as a function of the energy dealignment 
e. However, as the dephasing rate is increased, the oscilla- 
tions are gradually washed out and they essentially vanish 
in the strong dephasing limit with 7 = 0.4X2 /h. Thus, in this 
case, dephasing of the coherent oscillations of electrons in- 
side the DQD seems to reduce the signatures of interactions 
intheFCS. 




-4-2 2 4 



e[Q] 

Fig. 7 Factorial current cumulants with detector-induced dephasing. 
The parameters are Fd = (1/3)1^ and 7 = (black), 0.1^2/^ (red), 
0.2Q/h (puiple), OAQ/h (blue). The curves are shifted for clarity. 



7 Conclusions 

We have discussed our recent proposal to detect interactions 
among electrons passing through a nano-scale conductor by 
measuring time-dependent high-order factorial cumulants. 
For non-interacting electrons in a two-terminal scattering 
problem, the full counting statistics is always generalized bi- 
nomial, the zeros of the generating function are real and neg- 
ative, and consequently the factorial cumulants do not oscil- 
late as functions of the observation time or any other system 
parameter In contrast, oscillating factorial cumulants must 
be due to interactions in the charge transport. In cases where 
the factorial cumulants oscillate, the zeros of the generating 
function have moved away from the real-axis and into the 
complex plane. As we have shown, the motion of the dom- 
inant zeros of the generating function can be deduced from 
the oscillations of the high order factorial cumulants. 

Here, we have illustrated these ideas with a system con- 
sisting of two coherently coupled quantum dots attached to 
voltage-biased electronic leads. The dynamics of the DQD 
was described using a Markovian generalized master equa- 
tion which allowed us to treat strong coupling to the leads 
together with the coherent evolution of electrons inside the 
DQD. Interestingly, we found that even in cases where the 
Fano factor of the transport is sub-Poissonian, the high-order 
factorial cumulants still enable us to detect interactions among 
the charges passing through the DQD. Finally, we discussed 
the influence of detector-induced dephasing on the FCS and 
found that, for this model, such dephasing processes may 
reduce the oscillations of the high order factorial cumulants. 

Our work leaves several open questions for future re- 
search. It is still not clear exactly under what conditions in- 
teractions cause oscillations of the factorial cumulants. This 
will require further careful investigations of the singularities 
of generating functions in FCS, for example as in the recent 
work on singularities in FCS for molecular junctions IHTI . 
The answer to this question may moreover come from future 
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measurements of oscillating factorial cumulants in interact- 
ing nano-scale conductors. Finally, a new and promising di- 
rection of research combines the zeros of generating func- 
tions and high order statistics with dynamical phase transi- 
tions in stochastic many -body systems 14811491 . 

Acknowledgements We acknowledge many useful discussions with 
Markus Biittiker. The work was supported by the Swiss NSF. 



References 

1. L.S. Levitov, G.B. Lesovik, JETP Lett. 58, 230 (1993) 

2. Ya.M. Blanter, M. Buttiker, Phys. Rep. 336, 1 (2000) 

3. Yu.V. Nazarov (ed.), Quantum Noise in Mesoscopic Physics 
(Kluwer, Dordrecht, 2003) 

4. B. Reulet, J. Senzier, D.E. Prober, Phys Rev. Lett. 91, 196601 
(2003) 

5. W. Lu, Z. Ji, L. Pfeiffer, K.W. West, A.J. Rimberg, Nature 423, 
422 (2003) 

6. Yu. Bomze, G. Gershon, D. Shovkun, L.S. Levitov, M. Reznikov, 
Phys. Rev. Lett. 95, 176601 (2005) 

7. J. Bylander, T. Duty, P Delsing, Nature 434, 361 (2005) 

8. T. Fujisawa, T. Hayashi, R. Tomita, Y. Hirayama, Science 
312(5780), 1634 (2006) 

9. S. Gustavsson, R. Leturcq, B. Simovic, R. Schleser, T. Ihn, 
P. Studerus, K. Ensslin, D.C. Driscoll, A.C. Gossard, Phys. Rev. 
Lett. 96, 076605 (2006) 

10. C. Fricke, E Hohls, W. Wegscheider, R.J. Haug, Phys. Rev. B 76, 
155307 (2007) 

11. E.V. Sukhorukov, A. Jordan, S. Gustavsson, R. Leturcq, T. Ihn, 
K. Ensslin, Nature Phys. 3, 243 (2007) 

12. A. Timofeev, M. Meschke, J. Peltonen, T. Heikkila, J. Pekola, 
Phys Rev. Lett. 98, 207001 (2007) 

13. G. Gershon, Y. Bomze, E.V. Sukhorukov, M. Reznikov, Phys. Rev. 
Lett. 101, 016803 (2008) 

14. S. Gustavsson, R. Leturcq, M. Studer, 1. Shorubalko, T. Ihn, K. En- 
sslin, D.C. Driscoll, A.C. Gossard, Surf. Sci. Rep. 64, 191 (2009) 

15. J. GabeUi, B. Reulet, Phys. Rev. B 80, 161203(R) (2009) 

16. C. Flindt, C. Fricke, F. Hohls, T. Novotny, K. Netocny, T. Brandes, 
R.J. Haug, Proc. Natl. Acad. Sci. USA 106, 10116 (2009) 

17. C. Fricke, F Hohls, C. Flindt, R.J. Haug, Physica E 42, 848 (2010) 

18. C. Fricke, F. Hohls, N. Sethubalasubramanian, L. Fricke, R.J. 
Haug, Appl. Phys. Lett. 96, 202103 (2010) 

19. N. Ubbelohde, C. Fricke, C. Flindt, F Hohls, R.J. Haug, Nat. 
Commun. 3, 612 (2012) 

20. D. Kambly, C. Flindt, M. Buttiker, Phys. Rev. B 83, 075432 (201 1) 

21. C.W.J. Beenakker, H. Schomerus, Phys. Rev. Lett. 93, 096801 
(2004) 

22. H.F Song, C. Flindt, S. Rachel, 1. Klich, K. Le Hur, Phys. Rev. B 
83, 161408(R) (2011) 

23. H.F Song, S. Rachel, C. Flindt, 1. Klich, N. Laflorencie, 
K. Le Hur, Phys. Rev. B 85, 035409 (2012) 

24. Y. Komijani, T. Choi, F. Nichele, K. Ensslin, T. Ihn, D. Reuter, 
A.D. Wieck, arXiv: 1302.2200 (2013) 

25. A.G. Abanov, D.A. Ivanov, Phys. Rev. Lett. 100, 086602 (2008) 

26. A.G. Abanov, D.A. Ivanov, Phys. Rev. B 79, 205315 (2009) 

27. D. Kambly, D.A. Ivanov, Phys. Rev. B 80, 193306 (2009) 

28. D. Kambly, C. Flindt, M. Buttiker, J. Phys.: Conf. Ser 400, 
042026 (2012) 

29. G. KieBlich, P Samuelsson, A. Wacker, E. Scholl, Phys. Rev. B 
73, 033312(2006) 

30. A. Braggio, C. Flindt, T. Novotny, J. Stat. Mech. p. P01048 (2009) 

31. F Hassler, M.V. Suslov, G.M. Graf, M.V. Lebedev, G.B. Lesovik, 
G. Blatter, Phys. Rev. B 78, 165330 (2008) 



32. M. Kendall, A. Stuart, The Advanced Theory of Statistics (Charles 
Griffin & Company Limited, 1977) 

33. N.L. Johnson, A.W. Kemp, S. Kotz, Univariate Discrete Distribu- 
tions, 3rd edn. (Wiley, 2005) 

34. G. KieBlich, E. Scholl, T. Brandes, F Hohls, R.J. Haug, Phys. Rev. 
Lett. 99, 206602 (2007) 

35. R.B. Dingle, Asymptotic expansions: their derivation and inter- 
pretation (Dover, Washington, 1973) 

36. M.V. Berry, Proc. R. Soc. A 461, 1735 (2005) 

37. M.B. Plenio, PL. Knight, Rev. Mod. Phys. 70, 101 (1998) 

38. Yu. Makhlin, G. Schon, A. Shnirman, Rev. Mod. Phys. 73, 357 
(2001) 

39. S.A. Gurvitz, Ya.S. Prager, Phys. Rev. B 53, 15932 (1996) 

40. S.A. Gurvitz, Phys. Rev. B 56, 15215 (1997) 

41. C. Flindt, T. Novotny, A.P Jauho, Europhys. Lett. 69(3), 475 
(2005) 

42. C. Flindt, T. Novotny, A. Braggio, M. Sassetti, A.P. Jauho, Phys. 
Rev. Lett. 100, 150601 (2008) 

43. C. Flindt, T. Novotny, A. Braggio, A.P Jauho, Phys. Rev. B 82, 
155407 (2010) 

44. P Barthold, R Hohls, N. Maire, K. Pierz, R.J. Haug, Phys. Rev. 
Lett. 96, 246804 (2006) 

45. R. Sanchez, S. Kohler, P Hanggi, G. Platero, Phys. Rev B 77, 
035409 (2008) 

46. J. Zamastil, E Vinette, J. Phys. A: Math. Gen. 38, 4009 (2005) 

47. Y. Utsumi, O. Entin-Wohlman, A. Ueda, A. Aharony, Phys. Rev. 
B87, 115407 (2013) 

48. J.P Gan-ahan, 1. Lesanovsky, Phys. Rev. Lett. 104, 160601 (2010) 

49. C. Flindt, J.P Garrahan, Phys. Rev. Lett. 110, 050601 (2013) 



